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Abstract 

We present here a variational method for maximizing the bandgap in a one-dimensional 
system where the potential is subject to given constraints. Two specihc examples are studied 
in detail. In the first, we show that if the potential is constrained to lie between two values, the 
largest bandgap is obtained by a mixture of the highest and lowest potential - an exact result 
valid in any dimension. The second example fixes the first and second moments of the potential 
and seeks to extremize the bandgap. An exact result is obtained. Finally, we indicate how our 
techniques may be applied to photonic bandgaps. 


1 Introduction 

Over the past decade, the design and study of photonic bandgap materials have attracted a lot 
of attention Q. The design process often involves extremization of a bandgap with respect to 
variations in the dielectric constant profile in the material, subject to certain constraints. Materials 
so designed are expected to find use in optical communication. The theoretical work in this field 
lies in two broad areas: exact results that prove existence and provide bounds and numerical 
algorithms that provide optimal designs j^. In this paper we present a variational argument that 
provides criteria that the optimal design must satisfy, illustrated by two examples. (The calculations 
presented here are for quantum-mechanical bandgaps; however, the results are applicable also to the 
scalar Helmholtz equation - relevant for photonic band structure calculations - as indicated later.) 
The first example leads us to a maximum contrast rule: if the potential is allowed to vary over a 
range of values, the optimal configuration is a combination of the largest and smallest values. The 
second example starts with a constraint on the second moment of the potential energy and leads 
to a sufficient (but not necessary) relationship between the wavefunction and the potential that - 
rather surprisingly - are satisfied by the solution we obtain. We verify the criteria and examine their 
consequences in one-dimension. Our exact results may be used as benchmarks in testing numerical 
algorithms and as tests for sharpness of proposed bounds. 

We begin with the Schrodinger equation [Q 

— -I-ui/; = ei/i. (!) 

Upon imposing appropriate boundary conditions on the wavefunction ip, this selects the energy 
eigenvalues e„ and eigenfunctions ipn- We assume v to be real, and thus we can choose ipn to be real. 
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The index n = 0,1, 2, ... orders the energy levels, so Cn+i > £»• In one dimension with tp real, this 
index n can be taken to be the number of nodes in the wavefunction, or equivalently the winding 
number in the ip — ip' phase plane. In this paper - especially when we perform explicit calculations 

- we shall henceforth assume one dimension unless otherwise stated. 

We now wish to consider the following sort of question: What potential v - subject to constraints 

- gives a maximum e„ for a particular state n? More elaborate questions will be considered shortly. 

What type of constraints on the potential are we interested in? Since adding a constant v{r) —> 
v(pf) +uo gives the trivial extension e„ ^ + uq, the first constraint on v should probably be to fix 

the average of u, so 

{v) = J dTv{r) = vi- ( 2 ) 

The next constraint should limit the overall deviations from the average, and two reasonable choices 
might be to fix: 

1. the second moment, so 



2. the maximum contrast, so 

'^max ^ u(T) ^ Vmin- (4) 

In general, we assume a set of J constraints of the form 

Gj [u] = J dr gj [u] = Vj, j = 1,2,..., J. (5) 

(The maximum contrast constraint can be cast into this form by an appropriate limiting 
however, as we shall see, there are more direct methods.) 

Now suppose we have solved the Schrodinger equation ( 0 )- We now make a small perturbation 
of the form 

v{f) v{f) + Sv{f). ( 6 ) 

We evaluate the effect of this perturbation on the energy eigenvalue, to first order: 

Sen « {Sv)n- (7) 

If we seek e„ to be an extremum then we require our original state to satisfy 

0 = {Sv)n = J dr Svipl (8) 

for any allowed variation 6v. As an example of constraints, let us fix the first and second moments 
of v; then we require 


Thus, (||) can be satisfied in general if 


or 


0 = 

J dr 6v 

(9) 

0 = 

J drvSv. 

(10) 

ipl = a + Pv, 

(11) 

V = a + bip'^. 

(12) 
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This allows us to close the equations by solving the non-linear Schrodinger equation 


- V'^'lpn + [a -I- = enV'n 


( 13 ) 


The solution of the Schrodinger problem with the maximum contrast constraint is even easier, 
since we see from (H) that to maximize e„ we must choose 


r Vjnax, where > (j)^ 
\ Vrrnn, where 


(14) 


The constant (j) is chosen by imposing the first moment constraint. 
In general, then, {6v)n = 0 is satisfied by choosing 




.7 

= E« 


dg^ 

' dv 


(15) 


This equation can be inverted for v, giving v = where a represents the set of J parameters 

aj. We then solve appropriate non-linear equations: 

- + h[lljl\a]tljn = entpn- (16) 


In this paper, we wish to maximize not the energy eigenvalues of the Schrodinger equation itself, 
but rather the energy bandgap between two levels, subject to constraints on the potential similar to 
previously. However, the constraints on the average of v is not really necessary, since we consider only 
bandgaps. (Note that the indices on tp no longer refer to their ordering. In fact, we usually require 
that 'ip 2 -, '01 refer to consecutive eigenstates tpn+iy '4’n^ and so for periodic boundary conditions in 
one dimension, these wavefunctions - the band edges - will have exactly the same number of nodes.) 
Thus, requiring the bandgap to be an extremum with respect to constrained variations of u, we now 
have the following three coupled equations: 


- ^ 7^02 + Vi’2 
-I- Vipi 

-0? 


£202 

£101 


,7 

E 




dv 


(17) 

(18) 

(19) 


Again, we can invert to hnd v, and rewrite the problem as two coupled non-linear Schrodinger 
equations 

- V^02 -f /l[0i - 0i|a]02 = £202 (20) 

-V^01-b/i[0i - 0i|a]0i = £101. (21) 

These equations are no longer quite so simple. We now consider our previous examples, in one 
dimension. 


2 Two examples 

2.1 Maximum contrast constraint 

Without loss of generality, we reformulate the maximum contrast constraint as vq > v{x) > 0. Then 
returning directly to (Q), we find 

( 5(£2 — £i) = f dx [02 — 0i](5u < 0 ( 22 ) 
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with no constraint on the average of v. Clearly the solution is 


1 0 , < pjlix). 


( 23 ) 


Furthermore, we fix v{x) to be periodic with period L, and designate ipi and 'ip 2 as each having one 
node per period (of v(x)). In other words, the two selected wavefunctions are at the band edges 
with Bloch wave vector k = ^. The sole adjustable parameter left in the potential is the width of 
the ‘barrier’ of height uq, which we call 2 A {0 < A< L/2). Thus, we have: 

, , f 0 : -A<x <A 
= | no : A<x<L-A 


extended periodically. Symmetry dictates that ipi has nodes at a; = mL and ip 2 , at x = (n + ^)L 
(or vice versa) where m, n are integers. Therefore, 


and 


ipiix) 


cos(/cia:) 
oi sinh[Ki(a; — -j)] 


-A< X < A 
A < X < L - A 


'lp2{x) 


sin(fc2x) 
02 cosh[/t 2 (a: — f)] 


-A<x < A 
A < X < L — A 


both extended periodically. Here k'^ + k'^ = vq and Oi are constants that ensure matching at the 
boundaries and may be eliminated to give: 


2/1 tan 2/1 



2 / 1 ^ coth 




(24) 


and 


tan 2/2 
2/2 



coth 




(25) 


where 2/i = kiA, a = L/2A and ij = yToL/2. The variables to be solved for are yi,y 2 and a. A 
third equation is provided by (|^) which implies that 'tpi{x) = i/’|(a;) at a: = A (and x = L — A): 


cot^ yi 


(«-l)2/i^ \ 

^-2/i" ) 


= tan^ 2/2 1 


(a-l)2/2^ \ 
^-2/2^ / 


cot 2/1 
2/1 

tan 2/2 
2/2 



(26) 


The solutions to these equations are displayed in the figures. Fig-0 is a simultaneous plot of pi{A), 
/ 92 (A) and (e 2 — ei)T^/4 against a at ly = 5. Consistent with expectation, the energy difference is an 
extremum (maximum, in this case) when pi (A) = P 2 {A\ at a = 3.136; Fig.0is a simultaneous plot 
of ipi, '02 and v{x) for this set of parameter values. Fig.g shows rf{= uoT^/4) against 1 /q;(= 2A/L). 
The breaks in the plot are a consequence of poor convergence of the numerical routine used as K 2 
and then ki switch from imaginary to real with increasing rj. Fig.0 shows eiL^/4, e2T^/4 and (for 
reference) VqL^/A: against VqL'^ jA. Finally, Fig.^ shows the energy difference (times T^/4) in this case 
as well as for the potential v{x) = Vg cos^ ^ for the same Bloch states, for comparison. (The latter 
has analytical eigenfunctions in terms of Mathieu functions.) As might be expected, the difference 
is larger with our optimized choice of potential. 
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Figure 2: pi, p 2 and v {vo = 1-5) (dot-dash, dash and solid lines, respectively) against 2a:/L. 
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Figure 5: Size of bandgap (times i^/4) against for the optimized periodic square-well (solid 

line) and a sinusoidal potential with the same period and contrast (dashed line). 


2.2 Second moment constraint 


The equations we are to solve are 


■02 + '(''02 

= £202 

(27) 

01 + VTpl 

= £101 

(28) 

01 = av. 


(29) 


We proceed with a roundabout strategy to solve these coupled equations by first solving the pair of 
coupled equations: 


— 0''' + vip = eip 

(30) 

0^ = f3v. 

(31) 

This is equivalent to the single equation 

0 " - 0^/0 -1- e0 = 0, 

(32) 


which is is one dimension. Suppose we find two such solutions with the same u, yet independent. 
Then we also have a solution to (^), with a = (32—Pi- Clearly this gives a solution, but is it possible? 
The answer, as we demonstrate explicitly below, is yes. 

We anticipate a periodic potential (with period L) and further impose the condition that the two 
energy eigenstates we seek are at the band edges, with period 2L for both. (The eigenfunctions may, 
then, be chosen to be real.) The solutions to (^^ have remarkable properties , one of which is 
exploited here. The equation may be solved in terms of the (periodic) Jacobian elliptic functions sn 
or cn. (The third possibility, dn, turns out not to be appropriate.) We may now specifically require 
that the period of the potential be L. We finally obtain: 


i’2{x) 

(e2 — ei)L^ 
4 


k , 2Kx . 

, sn( ——) 

^L{1 - E/K) L 

k , 2Ka:. 

, cn( ——) 

^L{1 - E/K) L 

{kKf 


(33) 
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Figure 6: Size of bandgap (times ?^/4) against uL?' for the optimized potential (solid line), a 
sinusoidal potential with the same period and second moment (dot-dashed line) and the periodic 
square-well optimized (to a different criterion) in the previous example (dashed line). 


where E and K are the complete elliptic integrals of modulus k. Note that v is linearly related to 
both and '(/'i- Defining cr^ = D dx (w^ — {v)'^), we have 


aL^ 

It is, perhaps, useful to think of the whole system as being parametrized by k. The suitably scaled 
energy difference and strength of the potential (as a) may be easily obtained from equations (^) 
and (1^), respectively. Fig| shows (e 2 — ei)T^/4 against crL^. For comparison, the same quantities 
for a sinusoidal potential and the optimized potential in the previous example are also plotted. As 
expected, our optimized potential gives the largest energy difference. 


= 8K" 


2(1 -k fc2)(l - E/K) - P - 3(1 - E/K)2 


(34) 


3 Conclusion 

We have discussed strategies for bandgap extremization subject to two very different, but realistic, 
constraints on the potential energy function. Our conclusions, in both cases, are illustrated by 
one-dimensional examples. It is worthwhile pointing out that the variational argument in the first 
example applies explicitly in any dimension. The second relies on the serendipitous existence of two 
solutions to the equations ( p^ in one dimension. Whether this stroke of luck can be replicated in 
higher dimensions is, to our knowledge, an open question. 

We acknowledge useful discussions with Graeme Milton and David Dobson. This work was 
supported in part by a grant from the U. S. National Science Foundation. 
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